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Abstract : This paper focuses on the analysis of hydrodynamic 
loads on fixed offshore structures (horizontal cylinder) that are 
operating in shallow water and are often subjected to huge 
wave loading. For the purpose of this study, linear (Airy) wave 
theory was adopted together with the application of (21 ) in the 
load computation. The loads for six different sea states were 
computed using spread sheet for the following values of time t 
= 0, T/4, T/2. 
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1. Introduction 

Hydrodynamic wave loading on fixed offshore structures has 
been an issue of concern to the offshore oil and gas industry. The 
analysis, design and construction of offshore structures is 
arguably one of the most demanding sets of tasks faced by the 
engineering profession. Over and above the usual conditions and 
situations met by land-based structures, offshore structures have 
the added complication of being placed in an ocean environment 
where hydrodynamic interaction effects and dynamic response 
become major considerations in their design. In general, wave 
and current can be found together in different forms in the 
ocean. The existence of waves and currents and their interaction 
play a significant role in most ocean dynamic processes and are 
important for ocean engineers (Haritos, 2007) 

1.1 Hydrodynamics 

Hydrodynamics is concerned with the study of water in motion. 
In the context of an offshore environment, the water of concern 
is the ocean. 

Its motion, (the kinematics of the water particles) stems from a 
number of sources including slowly varying currents fromthe 
effects of the tides and from local thermal influences and 
oscillatory motion from wave activity that is normally wind- 
generated (Haritos, 2007). 
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2.1 Wave Theories 

All wave theories obey some form of wave equation in which 
the dependent variable depends on physical phenomena and 
boundary conditions (Al-Salehy, 2002). In general, the wave 
equation and the boundary conditions may be linear and non 
linear. 

2.2.1 Airy Waw Theory 

The surface elevation of an Airy wave amplitude £ a , at any 
instance of time t and horizontal position x in the direction of 
travel of the wave, is denoted by r|(x,t) and is given by: 

n(x,t)=^ a cos(kx-cot) (1) 

where wave number k = 27i/L in which L represents the 
wavelength (see Fig. 1) and circular frequency go = 27i/T in 
which T represents the period of the wave. The celerity, or 
speed, of the wave C is given by L/T or co/k, and the crest to 
trough wave height, His given by 2Q. 
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Figure 1 : definition diagram for an Airy wave (Haritos, 2007). 



2.0 Methodology 

The jacket structure used for this study is a HD accommodation 
platform to be operated in shallow water and is similar to all 
fixed jacket offshore structures. The part of the structure under 
water was discretized in to (264) Beam elements. The water 
depth for the HD field is approximately 25.3m. The loads were 
computed using spread sheet. 

Table 1: most probable wave heights and time periods for 
different sea states (Area 59) west Africa. 



the along wave u(x, t) and vertical v(x, t) water particle 
velocities in an Airy wave at position z measured fromthe Mean 
Water Level (MWL) in depth water h are given by: 

o)C, cosh l$k(z+h)] 

u(x,t)=— cosL'Pex - cot) (2) 

(jl)^ sinshiffi/c(z+/i)] 

v(x,t)= — sinL'Pcx - cot) (3) 

The dispersion relationship relates wave number k to circular 
frequency co (as these are not independent), via: 
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oo = gktanh (kh) 



(4) 



where g is the acceleration due to gravity (9.8 m/s ). The along 
wave acceleration it(x, t) is given by the time derivative of (2) 
as: 

a) 2 C cosh l>\k(z+h)] , . 

u(x,t)= ; sinUcx — cot) (5a) 



sinh ['(pfc/i ) 

while the vertical velocity x>(x, t) is given by the time derivative 
of (3) as: 

o) 2 ^ sinh['p(z+/i)] 

v(x,t)=- — coslfc - cot) (5b) 

v 7 sinh [fM) v y v ' 

It should be noted here that wave amplitude, a = c] a , is 
considered small (in fact negligible) in comparison to water 
depth h in the derivation of Airy wave theory. 
For deep water conditions, kh > tc, (2) to (5) can be 
approximated to: 

u(x,t)= co^e kz cosOcx — cot) (6) v(x,t)= a)C^e kz sinOcx — 

cot (7) 



co 2 =gk 



x>(x,t)= co 2 ^ e KZ sinl^lfex - cot) 



,kz 



(8) 
(9) 



This would imply that the elliptical orbits of the water particles 
associated with the general Airy wave description in (2) and. (3), 
would reduce to circular orbits in deep water conditions as 
implied by (6) and (7). 

2.2.2 Stoke 's Second Order Wave Theory 

Stokes employed perturbation techniques to solve the wave 
boundary value problem and developed a theory for finite 
amplitude wave that he carried to the second order. In this 
theory, all the wave characteristics (velocity potential, celerity, 
surface profile, particle kinetics... .e.tc) are formulated in terms 
of a power series in successively higher orders of the wave 
steepness (H/L). 

A condition of this theory is that (H/d) should be small so that 
the theory is applicable only in deep water and a portion of the 
immediate depth range. 

For engineering applications, the second -order and possibly the 
fifth-order theories are the most commonly used (Sorensen, 
2006). 

Stoke's wave expansion method is formally valid under the 
conditions (Iraninejad, 1988): 

H/d « (kd)2 for kd < 1 and H/L « 1 . 

Stoke's wave theory is considered most nearly valid in water 
where the relative depth (D/L) is greater than about (1/10) 
(Patal, 1989). 

Stoke's theory would be adequate for describing water waves at 
any depth of water. In shallow water, the connective terms 



become relatively large, the series convergence is slow and 
erratic and a large number of terms are required to achieve a 
uniform accuracy (Muga, 2003). 
The fluid particle velocities are then given by: 



TlH COSh[/c(z+/l)] rr.i fi ,n\ 
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(9a) 
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3(nH) 2 sinh [2k(z+h)] 
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(9b) 



The fluid particle accelerations are then given by: 

_ n 2 H cosh ik(z+h)] . r 1 \ 

a Y = 2 — ; — - — sinUcx — cot) + 



sinh (kh) 
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T 2 L sinh 4 (kh) 



sin2 (kx — cot) (10a) 
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These velocities and accelerations in (9) and (10) are used in 
Moris on 's equation to calculate load vectors of hydro dynamic 
loading by using Stoke's wave theory after being transformed 
from global coordinates for each member of the offshore 
structure. 

2.2.3. Moris on's Equation 

The along wave or in-line force per unit length acting on the 
submerged section of a rigid vertical surface -piercing cylinder, F 
(z, t), from the interaction of the wave kinematics at position z 
from the MWL 

(see Fig. 2), is given by Morison's quation. 
This equation is originally developed to compute hydrodynamic 
forces acting on a cylinder at a right angle to the steady flow, 
and is given by: 



F(z, t) = p7i — C m a(z, t) 
+ 7PDQ V\V\ (z,t) 



(11) 



in this equation, it is assumed that the wave force is acting on the 
vertical distance (z, t) of the cylinder due to the velocity (v) and 
acceleration (a) of the water particles, where (p) is the density of 
water, (D) is the cylinder diameter, (C m ) and (C d ) are inertia and 
drag coefficients, F D and F : are drag force and inertia force 
(Madhujit and Sinha, 1988). 

F D = 7PD C d V\V\ (z,t) (12) 



Fi= pTi— C m a(z,t) 

4 



(13) 
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These coefficients are found to be dependent upon Reynold's 
number, Re, Keulegan Carpenter number, KC, and the p 
parameter, Viz; 



D r KC 



(14) 




d [«h»L Cwrfimttt S,nhin 




(i) 

Figure 2: Water Particle Velocities along Member I [20] 
Where u m = the maximum along wave water particle velocity. It 
is found that for KC < 10, inertia forces progressively dominate; 
for 10 < KC < 20 both inertia and drag force components are 
significant and for KC > 20, drag force progressively dominates 
(Haritos,2007). 

The method adopted here assumes that only the components of 
water particles and accelerations normal to the member produce 
(Qian and Wang, 1992). 

To formulate the hydrodynamic load vector F w , consider the 
single, bottom mounted cylindrical member as shown in Fig. 2. 
The forces are found by the well known semi -empirical 
Morison's formula (11). It also represents the load exerted on a 
vertical cylinder, assuming that the total force on an object in the 
wave is the sum of drag and inertia force components. This 
assumption (introduced by Moris on) takes the drag term as a 
function of velocity and the inertia force as a function of 

acceleration (Zienkiewicz et al, 1978: Dean and Dairy mple, 
1984: McCormick, 1973), so that: 
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F n = P n ~ C 
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This can be simplified to: 



D 

Fn ~ P ft ~ Cm V n 



+ -pDC d (v n ). \(v n )\ 



(15a) 



(15b) 



Where: 

F n = nodal hydrodynamic force normal to the cylinder, D = 
Outer diameter of cylinder, 

p = Sea water density. C d = Drag coefficient (= 1.05). v„ = 
water particle acceleration. 

C m = Inertia coefficient (= 1.2). v n = water particle velocity. w n 
= Structural velocity. u^= Structural acceleration. 



Equation (15b) neglects the non -linear terms of drag coefficient 
(Al-Jasim, 2000; Sarpakaya and Isaacson, 1981). water particle 
velocity and acceleration can be evaluated by potential velocity 
computed from wave theories; the absolute value of velocity is 
needed to preserve the sign variation of the force. 

2.3 Global and Local System 

The kinematics of cross -flow with resultant velocity (see Fig. 3) 
is; 



[U L ]= [ti[u g ] 



(16) 



w n =^JTTTw 7 (17) 

is determined using wave theory applied in the global system 
and then transferred to the local system using transformation 
matrix. 




Global System 




Lutal System 



Figure 3: global and local system 



Application of Morison's equation leads to: 

17 D 2 _ du 

Ft = p 71 L m 

4 m dt 



+ 7PD C d (w n ) . |(w n )| 



(18) 



The components of the forces in the local axis system then 
become: 



= ipDC d (w n )l(w n )| 



u 

LW J 



D 



+ P 71 ~ C m 



U 
Lw J 



(19) 



To get the local forces, we need to get the member 
transformation matrix along the y - axis as follows; 
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Figure 4: global and local coordinates for a horizontal cylinder 
(y-axis). 
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Table 2: Hydrodynamic Loads @ time t = 0, distance x = 



And also the member transformation matrix along the x- axis is 
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Figure 5: global and local coordinates for a horizontal cylinder 
( x-axis) 
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Table 2: Hydrodynamic Loads @ time 
for (3) sea states 



t = T/4 and x = 16m 



Therefore, the local forces are given as:- 
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(20) 



The local forces are then transferred in to global forces by the 
transpose matrix 
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3.0 Tables & Results 

The following tables present the magnitude of hydrodynamic 
forces for three different sea states at different time phase. These 
forces serve as input loads for design. 



5.0 Conclusion 

It can be concluded that: 

The magnitude of all hydrodynamic forces along the member is 
uniform at constant depth. Also, the magnitude of the 
hydrodynamic force in the x-component increases in a regular 
pattern with change in sea state while the magnitude of the 
hydrodynamic force in the z-component increases in an irregular 
pattern with change in sea states. 
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